Tracking the effects of interactions on spinons in gapless Heisenberg chains 
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We consider the effects of interactions on spinon excitations in Heisenberg spin- 1/2 chains. We 
compute the exact two-spinon part of the longitudinal structure factor of the infinite chain in zero 
field for all values of anisotropy in the gapless antiferromagnetic regime, via an exact algebraic 
approach. Our results allow us to quantitatively describe the behaviour of these fundamental ex- 
citations throughout the observable continuum, for cases ranging from free to fully coupled chains, 
thereby explicitly mapping the effects of 'turning on the interactions' in a strongly-correlated system. 



Interactions in one-dimensional (Id) systems are 
known to overwhelm constituent particles, leading to 
a collective quantum liquid state with low-energy exci- 
tations described by the theory of Tomonaga-Luttinger 
liquids^. While the 'universal' physics of Id systems is 
phenomenologically well understood^, it almost always 
remains impossible to precisely track the effects of 'turn- 
ing on the interactions' on the constituent particles, as 
one does for Fermi liquids^ (where bare fermions are adia- 
batically connected to Landau quasiparticles) . In this re- 
spect, our general understanding of Id systems can bene- 
fit from nonperturbative solutions of microscopic models, 
a fundamental example being the Heisenberg spin- 1/2 
anisotropic chaiiP^, whose Hamiltonian is 

N 

H = jY: {Sp^^, + + A^I^I^O . (1) 

This system is a Tomogana-Luttinger liquid for 
anisotropy values A in the range —1 < A < 1 (in 
zero field, with J > 0). Its fundamental excitations are 
spinon^: spin-1/2 fractionalized objects which can be 
viewed as domain walls dressed by quantum fluctuations. 

A way to probe the nature of excitations is to deter- 
mine how they carry observable correlations, an interest- 
ing example here being the longitudinal structure factor 

At A = 0, this can be written as a density correlator of 
Jordan- Wigner fermions. Only single particle- hole exci- 
tations contribute, the exact structure factor being pro- 
portional to their density of states. For A > 0, this 
picture breaks dowri^ due to nonperturbative effects. 

It is the purpose of this paper to track in detail the ef- 
fects of 'turning on' interactions on the spinon quasiparti- 
cles and their ability to carry correlations, throughout the 
gapless antiferromagnetic regime < A < 1. Systems in 
this regime can be realized and studied experimentally 
(for fixed anisotropy) in spin ladder compounds^ or (in 
principle for generic anisotropy) using optical latticei^. 
Focusing on zero temperature, we will compute the exact 



two-spinon contribution to ([2]) directly in the thermody- 
namic limit N using an adaptation of the 'vertex 
operator approach™ Our results provide a strict lower 
bound and (for practical purposes) an extremely accu- 
rate representation for the complete correlator of the in- 
finite system (more that 99% for anisotropics below 0.5) 
throughout the observable excitation continuum. They 
provide a robust benchmark for assessing the linesh apes 
obtained for finite systems directly from integrability^!^^^ 
or using varia nts of the density matrix renormalization 
group (DMRCP^ or quantum Monte Carlo (QMC]P, 
and confirming the threshold behaviour predicted using 
field theor y^^^^ , complementing it with exact prefactors. 

The vertex operator approach was originally developed 
for A > 1 where the Hamiltonian commutes with the 
action of the quantum group Uq{sl2). The representa- 
tion theory of this quantum group leads to explicit ex- 
pressions for states, physical operators and their matrix 
elements^^, providing building blocks for correlations in 
terms of contributions from intermediate states made 
of increasing numbers of pairs of spinons, S^^{k^uj) = 
E^=i^f2m)(^'^)- The calculation of (g was treated 
using the vertex opera tor approach at A = 1 for twcP^^^ 
and four spinon^^^^, the combination being shown to 
yield about 99% overall accuracy. The A > 1 regime 
was also considered^^ . The physically more interest- 
ing quantum critical gapless regime (0 < A < 1) re- 
mains however largely unexplored by these exact ther- 
modynamic methods. Our paper aims to fill this gap. 

Spinon excitations - The ground state of the gapless 
XX Z antiferromagnet supports spinon excitations^ with 
exact zero-field dispersion relation e{p) = vf \ sinp|, p G 

[— 7r,0], where the Fermi velocity is i'f(A) = ^ acos^ • 
Spinons always appear in pairs, so the simplest states 
which contribute to the structure factor are made of 2 
spinons. Parametrizing their momentum by pi and p2, 
momentum and energy conservation impose k = —pi —p2 , 
uj = e{pi) -\- e{p2). The two-spinon states thus form a 
continuum in k-u defined by lower and upper boundaries 

^2,z(^) = '^f| sin a;2,n(^) = 2i'F sin /c/2. (3) 

Matrix elements via vertex operator approach - The 
vertex operator approach is also applicable, albeit in- 



2 



directly , to the gapless region < A < 1. The 
strateg3f^^^' is to first generalize the problem to the 
completely anisotropic Heisenberg model XK^^c^j '^j+i + 

3 

JySj Sj_^i + JzSj^Sj_^i) in the so called principal regime 

\Jy\ ^ Jx ^ Jz^ for which matrix elements of local op- 
erators between the vacuum and excited states can be 
computed exactly using a variant of the vertex operator 
approach^^*^. These results can then be mapped to the 
disordered regime \Jz\ < Jy < Jx^^ before taking the 
Jx Jy limit to reconstruct the matrix elements for the 
gapless Hamiltonian (1) with < A < 1. In this way we 
find^^ the following exact expression for the two-spinon 
contribution to S^^{k^w): 



Q{0J2,u{k) - UJ)Q{UJ - UJ2,l{k)) 



x(l + l/0' 



g-/5(p(fc,c^)) 



(4) 



COS ■ 



in which ^ 

h{p) ^ 



acosA 



1, 6 is the Heaviside function, and 
dt sinh(^ + l)t cosh(2t) cos(4pt) - 1 



/o t sinh^t cosht sinh(2t) 

in which the parameter p is defined as 



cosh(7rp(/c, uo)) 



''2,1 



(k) 



2,1 



(k) 



(5) 



(6) 



Results - In FigjT] we give plots of the two-spinon 
part of the longitudinal structure factor Q for values 
of A interpolating between weak and strong coupling. 
A few striking things are worth mentioning concerning 
the influence of interactions on the two-spinon part of 
the correlations. Most noticeably, upper threshold diver- 
gence disappears immediately upon turning interactions 
on. The correlation weight also starts flowing around the 
edges of the continuum, mostly via the wings at /c 0, 27r 
(see e.g. the A = 0.2 plot), and thereafter starts accu- 
mulating at the antiferromagnetic point k = tt (see the 
A = 0.4 plot). The lower threshold divergence starts car- 
rying more weight from A :^ 0.5 onwards, and becomes 
increasingly sharp as one approaches the isotropic point. 

Within the two-spinon continuum, so away from the 
thresholds, two things can be noticed. First, the weight 
within the bulk of this continuum quickly changes shape 
as A is turned on: from a pure [oJ2,u{k) — form 
at A = 0, it becomes almost uniform in frequency for 
A 0.2; it then becomes a rapidly decreasing function 
of frequency for higher interactions. Turning interactions 
on therefore leads to a remarkable collapse of correlation 
weight from high to low energies. 

Sum rules - To quantify the importance of the two- 
spinon contribution to the full structure factor, we use 
two useful sum rules, namely the integrated intensity 



Jo ^Jo 



duo 
2^ 



S(k,uj) = 1/4, 



(7) 



and the f-sumrule (at fixed momentum)!^, 

Il'{k)= I —uS{k,uj) = -2X''{l-cosk) (8) 
Jo 27r 

where = {SjSj_^i) is the ground state expectation 
value of the in-plane exchange term. This can be ob- 
tained from the ground-state energy density ej^ and its 
derivative, namely = ^(1 — A^)eo, with 

-J{^ 



1 



27r 



sm 



TT 



dt 1 



tanh t 



tanh[(^ + l)t] 



(9) 

We provide the explicit values of the sum rule saturations 
coming from two-spinon contributions in Table [l] (for the 
f-sumrule, the saturation is the same at all momenta). 
The two-spinon states carry the totality of the correlation 
at A = 0, and this remains approximately true up to 
surprisingly large values of interactions A ~ 0.8, above 
which four, six, ... spinon states become noticeable. 
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0.6 


0.9778 


0.9743 


0.1 


0.9997 


0.9997 




0.7 


0.9637 


0.9578 


0.2 


0.9986 


9.9984 




0.8 


0.9406 


0.9314 


0.3 


0.9964 


9.9959 




0.9 


0.8980 


0.8844 


0.4 


0.9927 


0.9917 




0.99 


0.7918 


0.7748 


0.5 


0.9869 


0.9849 




0.999 


0.7494 


0.7331 



TABLE I: Sum rule saturations as a function of anisotropy: 
two-spinon contribution to the integrated intensity ([7| 
and first frequency moment /f^ (|8|. 

Threshold behaviour - The behaviour of the longitu- 
dinal structure factor in the vicinity of the excitation 
thresholds can be determined from the analytic expres- 
sions we have obtained, allowing us to make contact with 
and complement recent field theory prediction J^^ * ^^l 

a. The structure factor near the upper threshold. 
The upper threshold uo ^2,u{k) is approached by the 
limit p ^ as can be seen from (|6|. A careful evalua- 
tion shows that the integral ([5| then behaves according 



to /^(p) 



21np + 0(1). We thus have from (6) 



that the structure factor vanishes as a square root. 



{k,uj) 



/„(C)(sin-)-^/yc.2,„(fc) 



(10) 

in which /n(0 ^ momentum-independent function 
of anisotropy. The anisotropy-independent square-root 
cusp at the threshold (for < A < 1) confirms the field 
theory predictions^^, and at A ^ 1 matches the same 
limit known to apply for the XXX casd^. The prefac- 
tor we obtain here varies quickly with momentum, show- 
ing strong enhancement of the upper threshold singular- 
ity when taking the momentum towards the k = 0, 27r 
zone boundaries (as can be seen in FigjT] most clearly at 
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FIG. 1: Two-spinon part of the longitudinal structure factor of the infinite Heisenberg chain, for different values of the anisotropy 
parameter A. For A ^ 0, the correlation follows the density of states, and has a square root singularity at the upper threshold 
for all values of momenta. Increasing the anisotropy shifts the weight progressively towards the lower boundary. The lower 
boundary becomes increasingly sharp as the A ^ 1 limit is approached. 
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small anisotropies). For the A ^ limit (so ^ ^ 1), the 
cosh + cos I in the denominator of |4| vanishes when 
p ^ 0. Overall, in this case one rather obtains a square- 
root divergence, Sn^ik.uo) — > f (l)-^^ ^^ 

which follows the singularity of the density of states (the 
matrix elements are then energy independent). This dis- 
continuous in A threshold exponent behaviour is also 
consistent with field iheor^'^^ . We notice further that 
the momentum dependence of the prefactor is changed 
to a much weaker one than that at A 7^ 0. 

h. The structure factor near the lower threshold. 
The limit uj Cj02,i{k) is obtained via p ^ 00. Evalu- 
ating (5) yields I^{p) > -tt (^1 + |^ P + 0(1). The 
structure factor then obeys 



> 



= M0 



\smk\ i)(sin|; 



1(1- 



(11) 

where /^(^) is again a momentum-independent function 
of anisotropy. The singularity exponent reproduces an 
early conjecture'^ and field theory predict ion^^^^^ the 
momentum-dependent part of the prefactor shows an 
even more complicated behaviour than that of the up- 
per threshold, being enhanced (though differently) both 
at the zone boundaries /c = 0, 27r as well as near k = tt. 
As a final detail, the A ^ limit (so ^ ^ 1) yields the 



expected behaviour, S2^{k^uj) 



> 



0(1). 



Conclusions - In this paper, we have tracked how the 
spinon excitations in Heisenberg antiferromagnets con- 
tribute to the longitudinal spin structure factor ([2|, as 
a function of anisotropy {i.e. interaction). We have ob- 
tained the two-spinon part of this correlator exactly in 
the zero-field, infinite-size chain throughout the gapless 
antiferromagnetic regime, by exploiting the vertex op- 
erator approach to express states and correlators in a 
purely algebraic language. Our results provide an exact 
lower bound for and an extremely accurate description 
of the full correlator (as shown by sum rule saturations) 
throughout the observable excitation continuum {i.e. not 
only at low energies or near thresholds), provide a re- 
silient check for alternate methods and give a nonper- 
turbative derivation of the threshold exponents obtained 
from field theory while complementing these with exact 
prefactor s. The precise functional form we have obtained 
also allows us to determine the region of validity of the 
threshold behaviour; we will address this and further is- 
sues in future worlP^. 
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